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Abstract 

Background: In woody plants from temperate regions, adaptation to the local climate results in annual cycles 
of growth and dormancy, and optimal regulation of these cycles are critical for growth, long-term survival, and 
competitive success. In this study we have investigated the genetic background to growth phenology in a Salix 
pedigree by assessing genetic and phenotypic variation in growth cessation, leaf senescence and bud burst in 
different years and environments. A previously constructed linkage map using the same pedigree and anchored to 
the annotated genome of P. trichocarpa was improved in target regions and used for QTL analysis of the traits. The 
major aims in this study were to map QTLs for phenology traits in Salix, and to identify candidate genes in QTL hot 
spots through comparative mapping with the closely related Populus trichocarpa. 

Results: All traits varied significantly among genotypes and the broad-sense heritabilities ranged between 0.5 and 
0.9, with the highest for leaf senescence. In total across experiment and years, 80 QTLs were detected. For individual 
traits, the QTLs explained together from 21.5 to 56.5% of the variation. Generally each individual QTL explained a low 
amount of the variation but three QTLs explained above 15% of the variation with one QTL for leaf senescence 
explaining 34% of the variation. The majority of the QTLs were recurrently identified across traits, years and environments. 
Two hotspots were identified on linkage group (LG) II and X where narrow QTLs for all traits co-localized. 

Conclusions: This study provides the most detailed analysis of QTL detection for phenology in Salix conducted so far. 
Several hotspot regions were found where QTLs for different traits and QTLs for the same trait but identified during 
different years co-localised. Many QTLs co-localised with QTLs found in poplar for similar traits that could indicate 
common pathways for these traits in Salicaceae. This study is an important first step in identifying QTLs and candidate 
genes for phenology traits in Salix. 
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Background 

Adaptation of the annual cycles of growth and dormancy 
to the local climate is critical for survival and competitive 
success of woody plants. Such local adaptation has been 
described in several species, and it can be seen as clines 
in phenology traits, including the timing of leaf emergence, 
leaf senescence, and growth cessation [1]. This pattern 
reflects the trade-off between frost tolerance and en- 
hanced growth [2]. Cessation and initiation of growth 
determines the period of active stem elongation, mark 
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the shift between frost resistant and vulnerable phases, 
and their timing is critical for overall biomass production, 
fitness and the long-term survival of species [3,4]. After 
height growth cessation and until leaves senesce, decidu- 
ous shrubs and trees, including willows and poplars, con- 
tinue to be photosynthetically active and may accumulate 
considerable biomass [5-7]. In this phase, photosynthesis 
is critical to cold acclimation and survival in winter [8]. 
The timing of leaf senescence in autumn has a strong 
impact on nutrient retranslocation, reserve storage and 
the next early-season growth [5,9]. In addition, spring 
and autumn leaf phenology are evolutionary important 
traits for herbivore and pathogen resistance [10-13]. 

Decreasing photoperiod (day-length) is the main envir- 
onmental cue inducing growth cessation and bud set in 
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many perennial plants [14], including poplar [15]. This 
response to photoperiod is under strong genetic control 
[16-18] and is maintained when trees are moved between 
latitudes [19]. Temperature, alone or in combination with 
photoperiod, also induces growth cessation in some tree 
species [20-24]. Senescence and shedding of leaves are 
also influenced by photoperiod, often in interaction with 
temperature [25]. In poplars, leaf senescence is induced by 
shortening day-lengths but needs to be preceded by bud 
set. However, it is not clear to what extent senescence and 
bud set are under independent genetic control [26,27]. 

After growth cessation, dormancy is initiated, which is 
a prerequisite for the development of cold acclimation 
and freezing tolerance [28]. Dormancy release requires 
exposure to chilling temperatures [29]. Light does not 
seem to play a role in this process, but once dormancy is 
broken, bud burst and growth resumption are regulated 
by temperature, light and photoperiod. Growth resumption 
and bud burst depend on accumulation of temperature 
units over a specific threshold (thermal time) [29]. Temp- 
erature is the most important factor regulating bud burst 
in temperate woody plants [30], but photoperiod also 
plays a role in some populations and species [31,32]. 

Willows (Salix) belong together with poplars (Populus) 
to the plant family Salicaceae. Based on the fossil record 
the divergence of the two genera was dated to approxi- 
mately 45 mya [33,34]. Willows and poplars share many 
characteristics such as dioecy, rapid growth and seed de- 
velopment, and ease with which they can be vegetatively 
propagated. They typically have a haploid chromosome 
number of 19 and similar genomes sizes of approximately 
500 Mbp. Also there is strong syntheny and colinearity 
between willow and poplar genomes [35]. The Salix 
genus shows a remarkable phenotypic diversity ranging 
from small shrubs to large trees. Salix spp. have a global 
distribution in temperate and arctic regions and are 
adapted to a wide range of habitats [36]. Relatively high 
levels of genetic diversity [37] and the broad phenotypic 
diversity make them an excellent model system for study- 
ing evolutionary processes such as adaptation. Moreover, 
willows have generally rapid growth and high biomass 
yields and these characteristics together with ease of vege- 
tative propagation make them economically attractive as 
bioenergy crops. Willows have been increasingly used in 
the last decades for biomass production worldwide and 
Salix viminalis L. and S. schwerinii E. Wolf and their 
hybrids are some of the most commonly used willows 
in the breeding programs in Europe. These two species 
are dioecious and outcrossing and morphologically very 
similar. Both are multi-stemmed shrubs with long and 
slender leaves and are commonly found along rivers and 
in other wet areas. 

Growth cessation in Salix species is marked by the 
abscission of the shoot apex [25] and is controlled by 



photoperiod [38,39]. Large clonal variation in the time of 
leaf abscission has been observed in Salix, and delayed 
leaf abscission was shown to impair leaf nitrogen retran- 
slocation and to increase nitrogen losses [6]. Extensive 
clonal and species variation in timing of bud burst has 
been observed in willow, mainly determined by differences 
in thermal time requirement [6,40,41]. Moderate to high 
heritabilities have been reported for timing of bud burst 
and growth cessation in different S. viminalis families 
[40,42,43]. 

Phenology traits have a quantitative genetic background 
and thus QTL mapping is a powerful method to identify 
genomic regions controlling phenology traits. With a 
reference genome one can obtain information on the 
genomic content of the QTL regions. With the advent 
of high-throughput genotyping technologies and anno- 
tated reference genomes, genetic markers in evenly spaced 
genes throughout the genome can be developed and geno- 
typed for the purpose of constructing dense genome-wide 
linkage maps. 

Here we studied the phenology of growth, including 
timing of bud burst, timing of cessation of elongation 
growth and leaf abscission in willows both in controlled 
and field conditions during multiple years. Our two aims 
were to map QTLs associated with phenology in Salix, 
and to identify candidate genes in QTL hot spots through 
comparative mapping with Populus trichocarpa. This was 
achieved by using a dense linkage map anchored to the 
annotated genome of P. trichocarpa and by constructing 
denser maps in the QTL hot spots. 

Results 

Phenotypic variation in phenology 

In the Si pedigree, planted in an experimental field in 
Pustnas, south of Uppsala (59°48' N, 17°39'E, 25 m), the 
mean date for bud burst was 20 th of April (day of the 
year (DOY) 112, Figure 1), the date of apex abscission 
(growth cessation) was 24 th of September (DOY 268, 
Figure 1), and about 25% yellow leaves were left on 
the plants at the end of October (leaf senescence index 
LSI = 1.5; DOY 304, Figure 1). Year-to-year variation was 
significant for all traits (Table 1, Figure 1), and genotype 
ranking significantly changed across years (Table 1, sig- 
nificant genotype x year interaction) indicating different 
responses to seasonal variation between genotypes. In 
the indoor experiment, elongation growth ceased on 
average 2 weeks after progressive reduction of the 
photoperiod. After nine weeks of artificial winter in the 
indoor experiment, bud burst occurred after 4 weeks of 
increased temperature and day length, corresponding 
to 336 day degrees > 0°C (Figure lg, h). 

In the progeny, all traits varied significantly among 
genotypes (Table 1, Figure 1). The variation of genotypic 
means for bud burst was smaller in the field than indoors 
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Figure 1 Distributions of phenotypic mean values for the S n pedigree which consists of 463 Ft progeny from the cross between the 
diploid hybrid male 'Bjorn' [Salix viminalis L. x S. schwerinii E. Wolf) and the diploid female S. viminalis '78183'. Distributions are shown 
for the traits bud burst (a 7 b), growth cessation (c, d) and leaf senescence (e, f) in field conditions during several years and indoors (g, h). 



(11 vs. 17 days, Figure lb, h). A large part of the pheno- 
typic variation for all the traits and in all experiments 
was due to genetic factors as shown by relatively high 
broad-sense heritabilities (0.5-0.9) (Table 2). The highest 



broad-sense heritability was found for leaf senescence. 
There were strong block effects (Table 1), but all geno- 
types were similarly affected in the field (Table 1, no 
genotype x block interaction) and indoors (among-blocks 
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Table 1 ANOVA for year and block x genotype interaction on phenology traits in the field 



Bud burst (BBF) 



df,/df 2 



Growth cessation (GCF) 
F df,/df 2 



Leaf senescence (LSF) 
F df,/df 2 



Genotype effect 


2.08*** 


309/309 


3.10*** 


139/139 


5.36*** 


452/452 


Year effect 


6405.57*** 


1/309 


1281.29*** 


1/139 


908.32*** 


1/452 


Genotype x year 


1 .95*** 


309/3100 


1 73*** 


139/1400 


2.47*** 


452/4530 


Genotype effect 3 


4.39*** 


309/1545 


5.57*** 


139/695 


13.99*** 


452/2260 


Block effect 3 


34.60*** 


5/1545 


8.16*** 


5/695 


66.03*** 


5/2260 


Genotype x block 3 


0.1 2 ns 


1545/1860 


0.26 ns 


695/840 


0.47 ns 


2260/2718 



a ANOVA performed on pooled data from 2 years. ***, Significant at the 0.01% level, ns= non significant. 



correlation r=0.5-07, p<0.001). Therefore, QTL analyses 
were performed on unadjusted clonal means, which 
should be a good estimate of the average behaviour of 
a genotype. 

Bud burst and growth cessation showed a significant 
correlation only in 2008. Positive but weak correlations 
between the field and indoor were found both for bud 
burst and growth cessation. In the field, bud burst, growth 
cessation, and leaf senescence were each positively cor- 
related between years (Table 3, plots in Additional file 1: 
Figure SI). 

QTL mapping 

In total across experiments and years, 80 QTLs were 
detected (Table 4, Figures 2 and 3). QTLs were named 
by the trait (BB = bud burst, GC = growth cessation, 
LS = leaf senescence), environment (I = indoors or F = 
field), year of assessment, linkage group where it was 
located, and if multiple QTLs were mapped to the same 
group, by an ordering number. On average about 40% 
(21.5-56.5% depending on the trait) of the variation 
among genotypic means was explained by a general model 
including all the QTLs that remained significant after 
backward selection. Individual QTLs generally explained a 
small proportion of the variance of the genotypic means. 



Eleven QTLs contributed to 7% or more of the trait 
variation, 3 QTLs explained more than 15%, and 1 
QTL explained up to 34% of the variation (Table 4). 

QTLs for bud burst 

For bud burst, 9 QTLs were identified indoors and 26 in 
the field during two years of assessments. Each of the 
QTLs included in the final model explained 1.0 to 14.1% 
of the variation in genotypic means (Table 4). The final 
model including all QTLs explained from 27.6% to 
53.2% of the observed variance depending on environ- 
ment and year of assessment (Table 4). No significant 
QTL by QTL interaction was detected for bud burst. 

QTLs for growth cessation and leaf senescence 

For growth cessation, 3 QTLs were identified indoors 
and 18 in the field during two years of assessments. The 
contribution of each QTL to the phenotypic variance 
was low to moderate (1.2% - 16.7%) but all QTLs together 
explained from 21.5 to 49.1% of the trait variation 
(Table 4). For leaf senescence, 24 QTLs were identified 
across the two years of assessment. The final model (10 
QTLs) explained 56.5% of the trait variation in 2009 
and 43.8% in 2010. The majority of QTLs made a low 
contribution to the total variance (1.0 - 6.0%) with the 



Table 2 Traits description, clonal mean heritability, and ANOVA results from variance components analysis 



Trait 



Year 



Name 



Unit 



F (d fi/df2) Genotype effect 



F (d fi/df2) Block effect 



Growth chamber experiment 












Bud burst 


2008 


BBI08 


Days 


0.78 


4.49**(293/586) 


1 2.25**2/586 


Growth cessation 


2008 


GCI08 


Days 


0.65 


2.86**(293/586) 


1 86.3** (2 /586) 


Field experiment 














Bud burst 


2009 


BBF09 


DOY 


0.52 


2.1 0**( 461 /2060) 


4-1 0 nS (5/2060) 


Bud burst 


2010 


BBF10 


DOY 


0.81 


5.30**( 461 /2415) 


136.4 nS (s/2415) 


Growth cessation 


2008 


GCF08 


DOY 


0.77 


4.30**( 461 /2157) 


1 1 .64^(5/2157) 


Growth cessation 


2009 


GCF09 


DOY 


0.71 


3.47**( 461 /2099) 


21 .84**( 5 /2099) 


Leaf senescence 


2009 


LSF09 


Senescence % 


0.92 


12.91**( 461 /2415 


18.73**( 5 /2415) 


Leaf senescence 


2010 


LSF 10 


Senescence % 


0.83 


5.86**( 461 /2296) 


91 -77^(5/2296) 



H 2 ,clonal mean heritability, **p<0,01, ns non significant. 
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Table 3 Pearson correlation coefficients between all phenology traits across years and experiments 


Ljruy Ljr i u 


DDCnO 1 

Doruy 


DDr I U 






DDI 




+0.70 


-0.01 


-0.09 


+0.40 


+0.37 


+0.14 


+0.20 


LSF09 














P<0 001 


P=0 82 


P=0 057 


P<0 001 


P<0 001 


P=0 02 


p=0 001 




+0.01 


-0.02 


+0.36 


+0.35 


+0.12 


+0.14 


LSF10 
















P=0.76 


P=0.65 


P<0.001 


P<0.001 


P=0.049 


P=0.02 






+0.35 


-0.001 


+0.07 


+0.19 


+0.05 


BBF09 


















P<0 001 


p=0 98 


p=0 87 


p=0 001 


P=0 39 








-0.04 


-0.04 


+0.30 


+0.11 


BBF10 




















P=0.37 


P=0.44 


P<0.001 


P=0.07 










+0.59 


+0.06 


+0.15 


GCF09 






















P<0.001 


P=0.29 


P=0.01 












+0.16 


+0.19 


GCF08 
























P=0.007 


P=0.001 














+0.02 


BBI 


























P=0.68 



^=462, 2 N=294 (see the Methods section for details). Correlations are based on genotypic means (3-6 individual plants per clone). 



exception of LSF09.II and LSF10.II, which explained 
34.2% and 19.3% respectively (Table 4). 

Comparison of QTL positions among years, environments 
and traits 

QTLs were considered to co-localize when their peak 
positions were less than 10 cM apart. Twenty-eight QTLs 
were identified only once while at 20 other positions at 
least 2 QTLs were mapped. Among the 20 regions where 
co-localizing QTLs were found, 11 included 2 QTLs for 
the same trait identified in 2 years or both environments, 
4 included 2 QTLs for different traits, and 5 included 3 
to 6 QTLs affecting different phenology traits (Table 4, 
Figures 2 and 3). 

For bud burst in the field, in five cases QTLs appeared 
more than once at the same map position comparing 
the two years (Figures 2 and 3). In contrast, for growth 
cessation in the field, most of the QTLs mapped to the 
same position both years (Figures 2 and 3). Six genomic 
regions were detected for leaf senescence where QTL 
appeared both years (Figures 2 and 3). 

In only one case were QTLs involved in bud burst found 
at the same map location when comparing indoors and 
the field conditions (Figure 2). One QTL for growth ces- 
sation was only identified indoors while another two 
co-localized with QTLs identified in the field (Figures 2). 

Improved mapping of two QTL hot spots for phenology 
on LG II and LG X 

The QTLs on LG II and X for bud burst, growth cessa- 
tion and leaf senescence were confirmed with the new 
denser linkage map. A comparison of the QTL mapping 
with the linkage map from Berlin et al [35] and the new 
maps is illustrated in Figures 4 and 5. The peak position 



of the QTLs was changed for all traits on LG II while 
only to some extent for the growth cessation traits on LG 
X. The 1.5 LOD intervals of the QTLs were considerable 
shortened with the new linkage maps (Figures 4 and 5). 
For the traits BBI and GCF09 the new analysis identified 
two QTLs on LG II for each trait instead of one as in the 
original analysis (Figure 4). 

Positional information of QTLs and identification of 
candidate genes 

We predicted the genomic interval for each QTL and 
summed all intervals for the traits, BBF09, BBF10, BBI, 
GCF08, GCF09, GCI, LSF09 and LSF10. This was possible 
by anchoring the SNP markers to the P. trichocarpa 
genome. The total genomic intervals varied between 
15.1 and 69.0 Mbp for the traits (interval for each QTL 
see Additional file 2: Table SI) and the total number of 
gene models ranged from 565 to 6,604. The largest 
number of gene models was found for BB10, however, 
39% is due to one QTL with a large 1.5 LOD interval 
on LG I. Since some QTLs covered the same genomic 
interval, some gene models appeared more than once, 
we therefore estimated the number of unique gene 
models for BB, GC and LSF to 9,633, 4,355 and 4,815 
respectively (Table 5). The total number of gene models 
in the intervals as well as the putative candidate genes 
for growth cessation are presented in Table 5. Several 
putative candidate genes were identified among these 
gene models (Figures 2 and 3; Additional file 3: Table S2). 
Candidate genes in the QTL intervals include photore- 
ceptors as well as several circadian clock genes and 
downstream components. Among photoreceptors both 
cryptochrome and phytochrome genes were identified 
(Figures 2 and 3; Additional file 3: Table S2). The phyto- 
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Table 4 Results of the QTL mapping procedure for each phenology trait 


QTL 


LG 


Marker 3 


P (cM) b 


SI (cM) c 


LOD 


a d 


PEV e 


Bud burst in indoor conditions (BBI) 














BB 


.lb 


lb 


l_21om_sa (L29gr.164f) 


150.9 


9.5 


7.86 


0.001 GW 


3.5 


BB 


.11.1 


llc?-2 


ll_30_sa_plll 


16.9 


23.5 


4.54 


0.005 CW 


7.5 


BB 


JI.2 


HcT-2 


ll_33_sa 


42.1 


17.5 


4.22 


0.007 CW 


EXCL 


BB 


JI.3 


ll<?-2 


ll_35_sa (L12b.203f) 


80.6 


12 


6.20 


0.006 GW 


7.1 


BB 


.VII 


VII 


Ro4b 


64.4 


48.5 


3.04 


0.032 CW 


EXCL 


BB 


.IX 


IX 


IX_1_sa 


56.6 


10 


3.60 


0.010 CW 


3.2 


BB 


.X.1 


X 


X-4 


83.7 


41 


3.55 


0.015 CW 


EXCL 


BB 


.X.2 


X 


X-21_sa (L3gr.227f) 


140.2 


19.2 


3.36 


0.023 CW 


2.8 


BB 


.A 


A 


XVI-13_sa 


17.0 


35 


3.95 


0.001 CW 


6.2 


BB 


total sum of explained variance 












30.3 


Bud burst in field conditions year 2009 (BBF09) 












BBF09.lb 


lb 


XVI_12_om_sa_pl 


99.5 


13 


6.65 


0.001 GW 


10.3 


BBF09.II 


HcT-2 


ll_24_sa 


59.9 


74.5 


2.78 


0.038 CW 


2.0 


BBF09.IV 


IV 


IV_20_sa (L2gr.226f) 


71.1 


87.9 


3.00 


0.028 CW 


EXCL 


BBF09.VI 


VI-1 


Vl_15_sa (L13y.371) 


48.6 


22.5 


5.24 


0.01 1 GW 


3.7 


BBF09.VII 


VII 


Vll-3b 


43.6 


70.5 


2.78 


0.043 CW 


EXCL 


BBF09.X 


X 


X_19_sa 


116.7 


32.5 


4.19 


0.003 CW 


1.4 


BBF09.XIV 


XIVcJ 


XIV_18_sa_plll 


85.9 


17 


5.89 


0.003 GW 


4.4 


BBF09.XVIII 


XVIII 


XVIII_6_sa 


14.6 


18.5 


3.25 


0.018 CW 


4.5 


BBF09.A 


A 


XVI_13_sa 


34.9 


18.5 


2.77 


0.001 CW 


EXCL 


BBF09.B 


B 


SB945 


27.7 


42.3 


2.28 


0.049 CW 


1.3 


BBF09 total sum of explained variance 












27.6 


Bud burst in field conditions year 2010 (BBF10) 












BBFIO.Ia 


la 


SB331 


170.3 


36.3 


4.60 


0.005 CW 


1.9 


BBF1 0.lb.1 


lb 


l_51om_sa 


107.2 


28 


3.78 


0.022 CW 


EXCL 


BBFIO.Ib.2 


lb 


l-32_sa 


215.5 


71.5 


3.63 


0.026 CW 


4.2 


BBF10.II 


HcT-2 


R_41_sa_pl 


20.9 


17 


13.32 


0.001 GW 


14.1 


BBF10.VI.1 


VI-1 


Vl-3d 


8.3 


15.5 


7.96 


0.005 GW 


7.3 


BBF10.VI.2 


VI-1 


SB496 


61.8 


13 


4.04 


0.005 CW 


EXCL 


BBF10.VII 


VII 


Ph18 


52.2 


28.5 


4.78 


0.002 CW 


1.0 


BBF10.VIII.1 


VIII-2 


Vlll_3_sa 


32.7 


29.3 


4.75 


0.001 CW 


5.6 


BBF10.VIII.2 


VIII-3 


Vlll_20_sa_pl 


41.4 


33 


6.72 


0.025 GW 


3.4 


BBF10.IX.1 


IX 


IX_5_sa 


36.9 


47.5 


3.79 


0.008 CW 


2.1 


BBF10.IX.2 


IX 


R_60_sa_pl 


78.2 


23.9 


3.73 


0.01 1 CW 


3.5 


BBF10.X 


X 


Ph2 


45.9 


25 


4.92 


0.002 CW 


3.2 


BBF10.XIV 


XIVc? 


XIV-4_sa 


82.3 


19 


3.38 


0.019 CW 


EXCL 


BBF10.XVII 


XVII 


ll_37_sa_pl 


123.8 


23.5 


3.19 


0.041 CW 


EXCL 


BBF1 0.XVIII 


XVIII 


XVIII_6_sa 


14.6 


28 


4.34 


0.002 CW 


3.0 


BBF10.A 


A 


XVI-13_sa 


35.0 


22 


3.05 


0.001 CW 


3.9 


BBF10 total sum of explained variance 












53.2 


Growth Cessation in indoor conditions (GO) 












GCI.Ib 


lb 


XVII_8om_sa 


182.8 


34 


3.93 


0.015 CW 


5.4 


GCI.II 


llc?-2 


ll_35_sa 


84.8 


14 


4.70 


0.006 CW 


4.6 
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Table 4 Results of the QTL mapping procedure for each phenology trait (Continued) 



GCI.V 




V-5_sa 


52.8 


24 


5.22 


0.028 GW 


11.5 


GCI total sum of explained variance 












21.5 


Growth cessation in field conditions year 2008 (GCF08) 












GCF08.la 


la 


R_29_sa 


42.0 


37.5 


3.38 


0.043 CW 


2.3 


GCF08.ID 


lb 


XVI_18_sa 


80.1 


12.5 


7.28 


0.001 CW 


6.1 


GCF08.II 


WS-2 


ll_35_sa 


84.8 


8 


10.54 


0.002 GW 


16.7 


GCF08.III 


III 


lll_8om_sa 


24.8 


36 


2.91 


0.045 CW 


1.7 


GCF08.V 


VcT 


V-5_sa 


52.8 


16 


5.93 


0.001 CW 


5.5 


GCF08.X 


X 


X-4 


90.4 


5.5 


16.31 


0.001 GW 


2.3 


GCF08.XIV 


XIVcJ 


XIV-4_sa 


82.3 


14 


4.42 


0.005 CW 


2.8 


GCF08.XVIII 


XVIII 


XVII_9om_sa_pl 


43.9 


47 


2.83 


0.043 CW 


3.5 


GCF08 total sum of explained variance 












40.9 


Growth cessation in field conditions year 2009 (GCF09) 












GCF09.ID 


lb 


XVI_18_sa 


80.2 


14.5 


4.75 


0.004 CW 


4.9 


GCF09.II.1 


HcM 


ll-lb_pl 


11.0 


26 


4.02 


0.001 CW 


7.0 


GCF09.II.2 


HcT-2 


ll_33_sa 


43.1 


17 


4.82 


0.004 CW 


2.4 


GCF09.II.3 


llc?-2 


ll_35_sa 


84.8 


13 


6.15 


0.001 CW 


6.2 


GCF09.III 


III 


lll-4_sa 


23.8 


29 


2.95 


0.047 CW 


3.2 


GCF09.V 


VcJ 


V_5_sa 


52.8 


55 


3.15 


0.032 CW 


1.2 


GCF09.IX 


IX 


IX_8_sa 


75.2 


22.9 


3.82 


0.007 CW 


4.2 


GCF09.X 


X 


X-4 


91.4 


5.5 


8.55 


0.002 GW 


10.7 


GCF09.XIV 


XIVcJ 


XIV-4_sa 


82.3 


7.5 


5.50 


0.001 CW 


1.4 


GCF09.XVIII 


XVIII 


XVIII_8_sa 


39.7 


51.5 


2.88 


0.042 CW 


2.1 


GCF09.lb x 


GCF09.XVIII 












5.8 


GCF total sum of explained variance 












49.1 


Leaf senescence in field conditions year 2009 (LSF09) 












LSF09.ID 


lb 


SB265 


219.5 


23 


4.34 


0.007 CW 


1.6 


LSF09.II 


HcT-2 


ll_35_sa 


84.8 


4 


32.11 


0.001 GW 


34.2 


LSF09.III 


III 


lll-14_sa 


87.1 


36 


3.57 


0.018 CW 


EXCL 


LSF09.IV 


IV 


Roll 


21.6 


38 


3.0 


0.030 CW 


EXCL 


LSF09.VI.1 


VI-1 


Vl-3d 


4.3 


10.5 


3.32 


0.023 CW 


1.2 


LSF09.VI.2 


VI-2 


Sa_con27_Popl 


8.4 


10.5 


8.11 


0.001 GW 


3.7 


LSF09.VII 


VII 


Vll_1_sa 


1.0 


6 


6.18 


0.001 GW 


4.3 


LSF09.VIII.1 


VIII-2 


Vlll_3_sa 


27.7 


23 


3.30 


0.005 CW 


1.6 


LSF09.VIII.2 


VIII-3 


Vlll_5_sa 


4.0 


15.5 


7.55 


0.001 GW 


5.6 


LSF09.X 


X 


X-4 


88.4 


6.5 


20.49 


0.001 GW 


2.2 


LSF09.XI 


XI 


R_52_sa 


4.5 


91 


3.14 


0.021 CW 


1.1 


LSF09.XIV 


xiv $ 


XIV-12_sa 


oU.U 


j4./ 


1 Al 


n n^n r\M 


1 n 
I .U 


LSF09 total sum of explained variance 












56.5 


Leaf senescence in field conditions year 2010 (LSF 2010) 












LSF10.II 


WS-2 


11-1 2_sa 


89.6 


5.5 


16.89 


0.001 GW 


19.3 


LSF10.III 


III 


lll-18_sa 


118.6 


15.5 


6.35 


0.003 GW 


6.0 


LSF10.IV 


IV 


Ph30 


95.3 


30.9 


2.88 


0.034 CW 


2.3 


LSF10.V 


VcT 


V_10_sa 


59.4 


28.5 


3.64 


0.005 CW 


5.0 


LSF10.VI 


VI-1 


Vl-3d 


4.3 


12.5 


2.94 


0.049 CW 


1.8 
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Table 4 Results of the QTL mapping procedure for each phenology trait (Continued) 



LSF10.VII 


VII 


Vll_1_sa 


1.0 


26 


2.73 


0.047 CW 


1.4 


LSF10.VIII 


VIII-3 


Vlll_5_sa 


4.0 


45.5 


2.91 


0.002 CW 


1.8 


LSF10.X 


X 


X-4 


88.4 


8.5 


10.49 


0.001 GW 


1.8 


LSF10.XII 


XII 


Xll-4d 


63.7 


45.5 


3.22 


0.016 CW 


1.2 


LSF10.XIV.1 


XIV? 


XIV-9_sa 


22.8 


20 


6.73 


0.001 GW 


3.2 


LSF10.XIV.2 


XIVcT 


XIV-14_sa 


135.8 


45.3 


4.60 


0.030 GW 


EXCL 


LSF10.B 


B 


Vl_22_sa 


12.0 


41.5 


4.51 


0.037 GW 


EXCL 


LSF10 total sum of explained variance 












43.8 



a The closest SNP marker to the QTL peak position, also used in multiple QTL mapping, and in mixed model and variance components analyses. When the closest 
marker to peak position was not a SNP, the AFLP marker is reported in brackets. b QTL position on the S q linkage map [35]. c Heuristic support interval determined 
as a 1.5-unit drop off on either side of the local LOD score peak. Probability for the null hypothesis of no QTL obtained after 5,000 permutations. Percentage of 
variance explained by the QTL or by the QTL x QTL interaction estimated by variance components analysis. GW Genome-wide, significant QTL. CW Chromosome-wide, 
suggestive QTL. EXCL QTL excluded from the complete model following a backward selection procedure at a test level of 0.05 (see Methods for details). 



chrome gene PHYB2 was located in a narrow range 
between two QTLs on LG X and should not be ruled 
out as a potential gene influencing growth control in 
willows. Gene models within QTLs also included several 
core circadian clock genes such as LATE HYPOCOTYL 
(LHY), PSEUDO-RESPONSE REGULATOR 7 (PRR7), LUX 
ARRHYTHMO (LUX), EARLY FLOWERING 3 (ELF3), 
ZEITLUPE (ZTL) and GIGANTEA (GI) (Figures 2 and 3; 
Additional file 3: Table S2) [44]. Interestingly the FT2 gene 
is located in the proximity of the fine mapped region 
on LG X (Figure 5). Noteworthy is that PtFT2 [45] and 
PttLHY [46] are involved in the control of growth cycle 
in Populus. 

Discussion 

The present study explores the genetic architecture of 
growth phenology in a pedigree between S. viminalis 
and S. schwerinii. Using a linkage map based on several 
hundreds of SNP markers, we identified QTLs for bud 
burst, growth cessation and leaf senescence in different 
years and environments. In the field, the QTLs explained 
together more than 40% of the variation in each trait. 
Several regions were identified where many QTLs co- 
localized for different traits or for the same trait across 
years. Since the SNPs were developed using the P. tri- 
chocarpa genome as a template, we obtained positional 
information for the Salix QTLs, projected them on the 
poplar genome, and identified the corresponding genomic 
intervals. The results suggest that some QTLs might be 
homologous to Populus QTLs. Moreover, in the projected 
QTL intervals we could identify putative candidate genes 
for the traits. 

Phenotypic variation 

All phenotypic traits varied among the progeny. A large 
part of the variation was imputable to genetic factors, 
which confirms moderate to high broad sense heritabilities 
for bud burst and growth cessation in Salix species 
[42,43] . Although heritability estimates are known to be 



environment and population specific, it is now well 
established that bud phenology is under strong genetic 
control in Salicaceae [16,18,47-51]. 

The behaviour of the progenies changed between con- 
trolled and field conditions. Growth cessation in short 
days indoors and apex abscission in the field displayed 
particularly weak correlations. One possible explanation 
to this discrepancy could be that the variation observed 
indoors only reflects a photoperiodic response, while the 
variation observed in the field could reflect effects of 
other environmental factors. In fact, short days alone can 
induce apical growth cessation both in seedlings and rooted 
cuttings of several Salix species, while in field conditions, 
apical growth cessation does not seem to be regulated by 
photoperiod, but other factors seem to be involved [38,52]. 
In S. viminalis, the coincidence of growth rate decline and 
tip senescence with development of low leaf water potential 
in summer suggests an effect by water stress [52]. 

The timing of phenological events changed across 
years, as shown by Weih [6] for other Salix clones. In 
addition, the ranking of the clones changed from one 
year to another, which indicates a plastic behaviour in the 
family. This might represent a differential response among 
the progenies to the seasonal differences in rainfall and 
temperature observed between years (Figure 6). Plasticity 
is crucial for a species to respond to the demands of a 
changing environment. An increasing number of studies 
indicate that temperature and stress factors may variably 
interact with photoperiod in controlling the timing of 
phenological events in woody species, including willow 
and poplar [20,53,54]. In S. viminalis, the effects of photo- 
period seem to be superimposed on those of water stress 
in controlling apical growth cessation [52]. A drier and 
warmer summer, i.e. higher temperature sum, higher max- 
imum temperatures and a greater number of dry days 
(Figure 6), might indeed be partly responsible for the 
markedly earlier growth cessation observed in 2009 in this 
study. Combined effects of photoperiod and temperature 
could be also responsible for the inter-annual variability 
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Figure 2 Location of spring and autumn phenology QTLs on the consensus genetic linkage map from the S-, pedigree on LG I to VII. 

The peak position of the local LOD (P in Table 4) is marked as a squared symbol and a support interval of 1.5-unit drop off on either side of the 
LOD peak is reported for each QTL QTLs are named as in Table 4. Candidate genes and their position in the corresponding Populus chromosome 
are indicated. Grey circles on LG II: Bud burst QTLs in field (BBF) and indoor (BBI) conditions from [42] (same mapping pedigree, tentative 
positions based on common AFLPs). QTLs for phenology traits in Populus (taken from literature) are indicated. Red boxes: Selected QTLs 
associated with bud set (several stages and sub-processes) in Populus nigra from Table six in [50]. White boxes: Bud set or bud flush QTLs from 
[48] at tentative positions, based on their distance from a SSR marker (Chr III) or phenology candidate genes (Chr VI) in the original publication 
(as in [18]) Green boxes: Bud set robust quantitative trait loci (QTL) regions in [18]. 
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(See figure on previous page.) 

Figure 3 Location of spring and autumn phenology QTLs on the consensus genetic linkage map from the S-i pedigree on LG VIII to XIX. 

The peak position of the local LOD (P in Table 4) is marked as a squared symbol and a support interval of 1 .5-unit drop off on either side of the LOD 
peak is reported for each QTL QTLs are named as in Table 4. Candidate genes and their position in the corresponding Populus chromosome are 
indicated. Grey circle on LG IX: Bud burst QTLs in field (BBF) conditions from [42] (same mapping pedigree, tentative positions based on common 
AFLPs). QTLs for phenology traits in Populus (taken from literature) are indicated. Red boxes: Selected QTLs associated with bud set (several stages and 
sub-processes) in Populus nigra from Table six in [50]. White boxes & Grey boxes: Bud set or bud flush QTLs from [48] at tentative positions, based on 
their distance from phenology candidate genes (Chr X) in the original publication (as in [18]) Green boxes: Bud set robust quantitative trait loci (QTL) 
regions in [18]. 



observed for leaf senescence in our willow pedigree. In 
Populus, yellowing of the leaves is initiated by a photo- 
periodic stimulus, but the progression of senescence is 
accelerated under low temperature [26]. Moreover in 
Populus, senescence seems to be faster in trees that 
have a late onset of senescence, independently on the 
effect of temperature. There is no data about regulation 
of seasonal leaf senescence in Salix. However, if we 



hypothesize a similar regulation as in poplar, the smaller 
total phenotypic variation (the full range of phenotypic 
variation in leaf abscission between individual plants was 
40 - 100% in 2009, while 10 - 100% in 2010; not shown) 
and the more advanced stage of senescence observed in 
2009 compared to 2010 could both be explained as an 
effect of lower temperatures (means, minima and maxima) 
in October 2009 (Figure 6a, b, c). 
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Figure 4 QTL hotspots on LG II for the mapped traits (BB = bud burst GC = growth cessation and LS = leaf senescence) and alignment 
to the poplar [P. trichocarpa) physical map. Homologous markers are connected with dotted lines. The positions of candidate genes involved 
in the photoperiodic pathway and circadian clock are shown on the poplar chromosomes. The regions with increased marker density are shown 
to the right of the arrows. 



Ghelardini et al. BMC Plant Biology 2014, 14:31 
http://www.biomedcentral.com/1471-2229/14/31 



Page 12 of 18 



PX 


s.,x 












s.,x 


PX 




_ 1 


0 


BBF10.X 

BBI.X.1 

GCF08.X 

GCF09.X 

LSF09.X 

LSF10.X 

BBF09.X 

BBI.X.2 




c 
c 


GCF08.X 
GCF09.X 
LSF09.X 
LSF10.X 


BBI.X.1 

GCF08.X 

GCF09.X 

LSF09.X 

LSF10.X 


80 cM 
- XJ6 

■ XJ39 •%*'*••.. 


16.6 Mbp 




TIC 
FKF1 _ 




Ph2 - 


"I 

- (.... 

■ *i 














- X_f36 


X_f6- 
XJ39- 

X_f35 
X f33- 




PFT1 - 

FT2 _ 
PRR3 _ 
T0E1 " 


- Ph2 
_ X-4 

- X-19 

. X-21 sa 

- 22.58 Mbp 


..•***' ...•**' x " 4 ■ 

•***" X_19_sa- 

X-21 sa. 

"\ 142.7 cm 11 


X-4 - 

Ro28_popl - 

X 1 7om sa . 
95 cM J 












: X j29 I:::::*—,,.;;; Xf29- 
X f20- 

.X fl7 

: j o f 2 i 8 6 - popl X - f17 " 

95 CM *a\ »lu 

18.1 Mbp 


: FT2 


Figure 5 QTL hotspots on LG X for the mapped traits (BB = bud burst GC = growth cessation and LS = leaf senescence) and 
alignment to the poplar [P. trichocarpa) physical map. Homologous markers are connected with dotted lines. The positions of candidate 
genes involved in the photoperiodic pathway and circadian clock are shown on the poplar chromosomes. The regions with increased marker 
density are shown to the right of the arrows. 





QTLs for phenology traits 

The phenology traits were typically quantitative. Several 
QTLs were found for each trait of which the majority 
explained less than 10% of the total variation. This was 
not surprising and seems to be a common feature of 
phenological traits among many species in the Salicaceae 
plant family [18,50,51]. Indeed, the regulation of phen- 
ology is quite complex and involves different pathways 
[14,55-57]. When estimating the total phenotypic variation 
of all QTLs for each trait we found particularly high 
values for leaf senescence, e.g. LSF10 explained as much 
as 56.5% of the total variation. The total phenotypic 
variation explained by the QTLs for the other traits 
was less strong but still substantial. 

Colocalization of QTLs 

At the level of resolution attained, there is considerably 
coincidence in map positions among a number of the 
QTLs. The most striking instances are found on LG II 
and X, where QTLs with narrow confidence intervals for 
BB, GC and LS consistently mapped across years and 
environments are located at similar positions. The QTLs 
on LG II also co-locate with bud burst QTLs found by 



Tsarouhas et al. [42] where they used a subset of the 
same willow pedigree and a different linkage map with 
some common markers (Figures 2 and 3). Fine mapping 
with additional markers on LG II and LG X further sup- 
ported that QTLs for all three traits might represent the 
same locus. Similarly, on LG V, QTLs for growth cessation 
and leaf senescence co-localized, on LG VI, QTLs for leaf 
senescence and bud burst co-localized and on LG XIV 
QTLs for bud burst and growth cessation were found at a 
similar genomic position. This suggests that pleiotropic 
effects of individual QTLs on several traits could result 
from shared components of the pathways controlling the 
different traits. This is perhaps not unexpected given that 
both bud burst and growth cessation are both controlled 
jointly by photoperiod and temperature. The time of leaf 
senescence has been shown to be an important trait for 
seasonal acclimation as well as for biomass production in 
willows [6]. It generally occurs well after growth cessation, 
but data from Populus tremula suggest that also the onset 
of leaf senescence is under photoperiodic control although 
it might respond to a different photoperiod than bud set 
[26]. The frequent co-localization of QTLs for growth ces- 
sation and leaf senescence suggests that genes controlling 



Table 5 Genomic intervals and number of gene models in QTL regions 



Trait 


No. of QTLs with positional information 


Total genomic interval of QTLs (Mbp) 


Range (Mbp) 


No. of gene models 


BBI 


8 


17.8 


0.3 -7.6 


1,841 


BBF09 


7 


30.0 


0.5-13.3 


2,676 


BBF10 


17 


69.0 


1 .4 -26.7 


6,604 


GCI 


5 


15.1 


1.4-11.8 


565 


GCF08 


10 


33.8 


0.04-11.8 


2,912 


GCF09 


11 


29.6 


0.04-11.8 


2,836 


LSF09 


4 


18.8 


0.04-4.3 


2,106 


LSF10 


5 


38.1 


0.04-11.8 


3,412 
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Figure 6 Meteorological variables at the field site in Ultuna (Sweden) during the three years of study and in relevant periods of time. 

Monthly minimum (a), mean (b), and maximum (c) temperatures; number of days with no rain (cumulated dry days, CMDD) since August 1st 
(d); cumulative sums of daily mean temperatures (CMT,°C) during August (e) and October (f); number of chill days (CD) with mean temperature 
below 5°C (g); cumulated mean temperatures above 0°C (day degrees>0°C) from March 1st to bud burst (h). 



Ghelardini et al. BMC Plant Biology 2014, 14:31 
http://www.biomedcentral.com/1471-2229/14/31 



Page 14 of 18 



growth cessation may to a large extent also affect the 
timing of leaf senescence. Furthermore, the common QTLs 
on some of the linkage groups suggest that these locations 
indeed contain genes with a central role in controlling 
seasonal growth in willows. 

Comparison of QTL regions between poplar and willows 

When comparing the willow QTLs with QTLs found in 
poplar many QTLs were found in similar genomic posi- 
tions, which suggests common mechanisms controlling 
these traits in willows and poplar. Rohde et al. [18] mapped 
six narrow QTLs for bud set in poplar to LG III, V, VI 
and two on VIII and XIII, and we find QTLs for growth 
cessation or leaf senescence at five of these positions 
(Figures 2 and 3). Fabbrini et al. [50] identified several 
QTLs for bud set in P. nigra of which those on LG I 
and IV overlap with QTLs for growth cessation and leaf 
senescence in willows. Frewen et al. [48] mapped three 
QTLs in poplar for bud set to LG III, VI and X of which 
those on III and VI was also identified in willows. There 
are also some major differences as for example the 
QTLs on LG II and IX in willows are missing in poplars. 
However the presence or absence of a QTL for a trait 
depends on several factors such as accuracy of the pheno- 
typing, environmental factors and perhaps most import- 
antly on whether or not the trait or gene is variable in the 
mapping population under study. In this study we have 
used one mapping population from a back cross between 
two species and how general these results are should be 
validated in another genetic background. 

Identification of positional candidate genes 

Comparative mapping with poplars represents a first 
step to further dissect the genetic basis of phenological 
variation in willows and ultimately identify genes or al- 
leles responsible for the trait variation that we observe. 
Our understanding of the physiology and biochemistry 
of the traits of interest very much define the success of 
this approach as it relies on previous identification of 
genes potentially involved in the control of the traits as 
well as on genome conservation between poplars and 
willows. We have previously shown that overall gene 
order is conserved between willows and poplars except 
for few large-scale chromosomal rearrangements [35], 
which justify anchoring the Salix QTLs to the physical 
map of poplar. We started by estimating the number of 
genes within the QTLs for each trait and found as 
many as 9,633 genes for bud burst and about half the 
number of genes for growth cessation and leaf senes- 
cence. Functionally characterized genes for phenology 
are still scarce in poplar apart from two examples. The 
poplar gene PtFT2 located on LG X that is known to 
be involved in growth cessation [45], was located close 
to the fine mapped region. LHY-genes have recently 



been shown to be involved in both growth cessation 
and budburst [46] and LHY1 was located close to the 
QTL hot spot region on LG II and LHY2 was located in 
a cluster of QTLs for bud burst and growth cessation 
on LG XIV. The fine mapped region on LG II contains 
a SVP homolog. In perennial species SVP genes have 
been shown to be involved in the growth cycle typically 
with high expression levels during bud dormancy [58-61]. 
The above-mentioned genes are strong candidate genes 
for growth cessation and bud burst that warrant further 
investigation in willows. 

Conclusion 

We identified substantial variation in all traits in the 
pedigree and all traits were associated with many QTLs 
that each explained less than 10% of the variation, a 
typical pattern of quantitative characters. In total, we 
identified 80 QTLs, of which some were clustered in 
hotspots where QTLs for the different traits co-localised. 
Two such hotspots on LG II and X were further investi- 
gated by the construction of denser linkage maps in these 
regions, an effort that greatly reduced the QTL intervals 
(and number of gene models). Some QTLs appear to co- 
localize with those found in poplars, which could indicate 
common pathways for these traits in Salicaceae. This 
study is an important first step in identifying QTLs and 
candidate genes for phenology traits in Salix but further 
work is needed e.g. to confirm the QTLs in other genetic 
backgrounds, further fine mapping and functional studies, 
to verify candidate genes. 

Methods 

Plant material, experimental design and phenotyping 

The Si pedigree consists of 463 F x progeny from the 
cross between the diploid hybrid male 'Bjorn (Salix vimi- 
nalis L. x S. schwerinii E. Wolf) and the diploid female S. 
viminalis 78183' originating from southern Sweden. The S. 
schwerinii parent (79069) of Bjorn originates from Siberia 
while the S. viminalis parent is an interspecific cross be- 
tween the male clone 78101 from Western Sweden and the 
female clone 78195 from southern Sweden. The parental 
clones of Si were selected based on variation in phenology 
traits [42]. The pedigree is conserved in an orchard near 
to Uppsala (59°49' N 17°40' E, central Sweden) where all 
the plant material used in this study was collected. 

Experiment 1 - Growth cessation and bud burst indoors 

In spring 2008, growth cessation was assessed in the S 2 
pedigree in a phytotron experiment under controlled day 
length and temperature conditions. The two parental 
genotypes and 294 randomly drawn genotypes were 
propagated by means of hardwood cuttings and planted 
in 1.1 litre pots filled with Weibulls 'Kron Mull' (organic 
matter 95%; pH 5-5-6.5; 180 g m" 3 N, 110 g m" 3 P, 
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195 g nT 3 K, 260 g m" 3 Mg, 100 g m" 3 S, 2000 g m" 3 Ca). 
All genotypes were assigned to each one of three walk-in 
growth chambers in a complete randomized block design, 
where each chamber represented a block with one rep- 
licate. Plants were grown for four weeks under 20°C 
constant temperature, 70% relative humidity and 20 h 
photoperiod (300 \imo\ PAR m" 2 s" 1 ). The day length 
was then reduced to 16 h for one week, to 14 h for 
another week, and then one hour per week down to 10 
h. After the first two weeks of growth, the plants were 
pruned and only the main shoot was preserved for the 
experiment. Plant height, i.e. the length of the stem 
from the emerging point on the wooden cutting to the 
tip of the apex, was measured once a week during week 
three and four, and three times per week thereafter. A 
sigmoid curve on the form Y=b/(l+a x e ( " k x x) ), where 
Y is plant height and X is the number of days since the 
first measurement, was fitted to the data of each plant 
in order to estimate the end of elongation growth, i.e. 
the day on which the estimated plant height reached 
95% of the final value. The date of growth cessation 
was expressed as the number of days since the start of 
the reduction of daylength. 

After the simulated autumn, the plants were subjected 
to an artificial winter at 8°C constant temperature and 9 
h photoperiod for 9 weeks. The plants were then cut 
back, leaving 5 cm stem. Bud burst was forced by keep- 
ing a constant temperature of 12°C for six weeks and 11 
h photoperiod. During this time the plants were checked 
for bud burst daily or every second day, depending on 
the speed of the process. Bud burst was defined as stage 
3 according to the phenological scale previously used by 
Weih [6]. Individual plants were recorded as flushing on 
the day when at least one bud reached stage 3. The date 
of bud burst was expressed for each plant as the number 
of days since the beginning of forcing. 

Experiment 2 - Bud burst, growth cessation and leaf 
senescence in the field 

In spring 2008, plants of 463 genotypes were planted in an 
experimental field in Pustnas, south of Uppsala (59°48 ' N, 
17°39'E, 25 m), at a spacing of 130 x 50 cm (i.e., about 
20,000 plants ha~ 1 ), according to a randomized complete 
block design comprising six blocks with one plant per 
genotype in each. Two border rows were planted around 
the experiment to reduce marginal effects. The plants 
were obtained by 5 cm hardwood cuttings rooted in 0.5 L 
peat pots with Weibulls Kron MulF as growing medium. 
The plants were grown for five weeks in a greenhouse and 
then transferred outside for hardening. Before planting, 
the site was appropriately prepared [62], including plough- 
ing, harrowing and repeated application of a systemic 
herbicide (Glyphomax, Dow AgroSciences, Indianapolis, 
IN). The plantation was irrigated in summer 2008 and 



weed controlled during the whole experimental period. 
The plants were cut back in winter 2009 and fertilized in 
spring 2009, 2010 with NPK (21-4-7) corresponding to 
80 kg N/ha and year. 

Growth cessation as defined by shoot apex abscission 
was scored using the highest shoot of each plant in 2008 
and 2009 from the end of August and once to twice 
per week depending on the rate of progression. Leaf 
senescence and abscission was visually estimated on 
October 31 2009 and 2010 according to the following 
leaf senescence index (LSI): 0 = no leaves left on the 
plant (100% abscission); 0.5 = less than 10% brownish 
leaves (~ 95% abscission); 1= 10 to 20% brownish 
leaves (~ 85% abscission); 1.5 = 20 to 30% brownish or 
yellow leaves (-75% abscission); 2 = 30 to 40% yellow 
leaves (-65% abscission); 2.5 = 40 to 50 yellow and 
green leaves (-55% abscission); 3 = 50 to 65 green 
leaves (-40% abscission); 3.5 = 65 to 80% green leaves 
(-30% abscission); 4=more than 80% green leaves 
(-10% abscission). Bud burst, defined as in Experiment 1, 
was assessed twice a week during April and May 2009 and 
2010. Date of bud burst and date of apex abscission were 
expressed as day of the year (DOY), i.e. number of days 
since January 1. Plants were cut back in January 2009. 
Therefore spring phenology was assessed on stumps in 
2009 and on one-year old shoots in 2010; and autumn 
phenology was assessed on one-year shoots in 2008 and 
2009, and on two-year old shoots in 2010. In order to 
characterize the weather conditions at the plantation 
site, temperature sum (cumulated mean temperatures 
CMT,°C), number of chill days with mean temperature 
below 5°C (CD<5°C), and the number of days without 
rain (cumulated dry days, CMDD) were calculated in 
relevant periods of time in all years of study from the 
records of a nearby meteorological station (Figure 6). 

Statistical analyses and QTL mapping 

The complete set of data from each experiment, which 
included all genotypes that had records for at least three 
ramets, was analysed with a mixed model ANOVA and 
variance components analysis to determine the effects of 
genotype (set as random factor) and block (set as fixed 
factor), and to estimate the broad-sense heritability, i.e. 
the ratio between genetic variance and total phenotypic 
variance. The genotype x block interactions in the field 
were tested using grouped data from two years. Subse- 
quent analyses were performed on unadjusted mean 
values among blocks in each experiment. Correlations 
among genotype means were calculated and plotted 
(Additional file 1: Figure SI) for all combinations of 
traits in each experiment and across experiments. The 
year effect and its interaction with genotype were tested 
in a separate ANOVA for all the traits assessed in the 
field. QTL analyses were performed with MapQTL @ 6.0 
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[63] and the linkage map previously developed for the 
Salix pedigree Si was used [35]. First, interval mapping 
[64] was applied using 1.0 cM steps across the genome 
to determine putative QTLs involved in the variation of 
each trait. Following Churchill and Doerge [65], the loga- 
rithm of the odds (LOD) threshold for QTL significance 
was empirically estimated from 5,000 permutations of 
phenotypic data. Two theoretical critical thresholds 
were considered for detection of a putative QTL: the 
first corresponding to genome wide error rate of 5% 
was used to define significant QTLs, and the second 
corresponding to a type I error of 5% at the chromosome 
level was used to define suggestive QTLs (suggestive link- 
age; [66,67]. Multiple QTL mapping (MQM) [68] was then 
performed on the same data: the nearest SNP marker to 
each putative QTL peak was used as a cofactor to control 
the genetic background while testing at a position in 
the genome. Only markers close to QTLs significant at 
the genome wide level were used as cofactors in MQM. 
When a cofactor was also a flanking marker of the 
tested region, it was automatically excluded from the 
model. The number of cofactors used varied between 1 
and 5. A 1.5-unit drop off on either side of the local 
LOD score peak was used to determine heuristic support 
intervals for significant QTLs. 

A backward selection procedure, at a test level of 0.05, 
was performed to determine whether a particular QTL 
could be dropped from the model resulted from the 
MQM analysis. For each trait, the model involved the 
genotype at the closest marker to the corresponding 
putative QTL. Pairwise epistatic interactions between 
all putative QTLs were tested, via the corresponding 
marker x marker interaction treated as random effects, 
by comparing the deviance of a model including all 
main effect QTLs and the specific epistatic interaction, 
with the deviance of a model including only the main 
effect QTLs and no epistatic interaction [69]. The deviance 
test (p < 0.001) was applied for testing interactions. 
The contribution of each significant QTL and epistatic 
interaction was then estimated by maximum-likelihood 
variance components analysis. Statistical analyses were 
performed using R [70] and SPSS v. 19. 

Improved mapping of two QTL hot spots for phenology 
on LG II and X 

Genomic DNA from 463 genotypes from the Si pedigree 
was extracted as in Berlin et al. (2010). Two chromosomal 
regions on LG II (75-90 cM) and LG X (84-93 cM ) 
containing QTLs for all traits except bud burst in the 
field 2010 were selected to increase the marker density 
and new primer pairs were designed in these regions 
(Additional file 4: Methods and Additional file 5: Table S3). 
Gene segments were amplified by PCR and sequenced in 
each parent as described in Berlin et al [35]. A total of 96 



SNPs were selected for genotyping, 59 from LG II region 
and 37 from LG X region. Linkage maps were constructed 
as in Berlin et al. [35], for details see Additional file 4: 
Methods. 

New QTL analyses for group II and X with the denser 
linkage maps were conducted. In the linkage map used 
for the QTL analysis all AFLP markers were removed 
since only 96 out of 463 individuals were genotyped with 
AFLPs and thus did not add much information to the 
analysis. When several markers were located at the same 
position only the most informative was kept. MQM 
mapping using MapQTL @ 6.0 [63] was conducted for 
comparison of QTLs between new and earlier linkage 
maps. The most informative marker close to the peak 
position of the QTL was used as cofactor in the MQM 
analysis. 

Comparative mapping and in-silico selection of 
candidate genes 

Physical coordinates of the QTLs were obtained from 
anchored markers in the P. trichocarpa genome assembly 
version 3 (http://www.phytozome.net/poplar). For each 
QTL interval, positions of sequences containing the SNP 
markers flanking 1.5-LOD on both sides of the LOD score 
peaks were determined by BLASTN searches. Once the 
regions were determined, gene models (predicted by 
the Gnomon gene prediction tool) were downloaded 
using BioMart that were subsequently annotated using 
the P. trichocarpa version 3 annotation information. 
Annotated gene models positioned in QTLs were searched 
for putative candidate genes. Since growth cessation 
and bud set are mainly controlled by photoperiod, 
genes in the photoperiod pathway and the circadian 
clock were considered as candidate genes for these 
traits (Additional file 3: Table S2). Genes controlling 
budburst and leaf senescence are less well known and 
candidate genes for those traits with functional charac- 
terization are scarce or absent. We therefore did not 
attempt to identify any such candidate genes. 

Additional files 



Additional file 1: Figure SI. Pairwise plots of phenology traits across 
years and environments based on mean values for each individual in the 
mapping population SI. Red line show the linear fit of points. 

Additional file 2: Table SI. Genomic regions for each QTL. 

Additional file 3: Table S2. Candidate genes within QTL regions for 
each trait. Positions of the genes are indicated in Figures 2 and 3. 

Additional file 4: Methods. Improved linkage map in two QTL hot 

spots for phenology on LG II and X. 

Additional file 5: Table S3. Markers in fine mapping areas. 



Competing interests 

The authors declare that they have no competing interests. 



Ghelardini et al. BMC Plant Biology 2014, 14:31 
http://www.biomedcentral.com/1471-2229/14/31 



Page 17 of 18 



Authors' contributions 

LG performed phenotypic measurements, conducted and interpreted the 
phenotypic and QTL analyses and wrote the manuscript. SB created the new 
linkage map, performed the genome analysis and wrote the manuscript. 
MW read and commented on the manuscript and provided funding for 
measurements in the field experiment. NG and UL were involved in 
phenotypic measurements, selected candidate genes, read and commented 
on the manuscript. ACRW conceived the study, performed phenotypic 
measurements, conducted QTL analysis with the new linkage map, wrote 
the manuscript and provided funding for the project. All authors read and 
approved the final manuscript. 



Acknowledgement 

We thank Nils-Erik Nordh, Richard Childs and Lorenzo Bonosi for conducting 
phenotypic mesurements in the field; Marika Pettersson and Andreas Fritz 
for phenotypic measurements indoors. We also thank Pascal Pucholt for 
constructing the pearl script for the primer design and Francois Besnier for 
suggestions on the QTL analyses. This project was supported by funds 
from the Swedish Energy Agency and the NL-faculty at the Swedish 
University of Agricultural Sciences. 

Author details 

department of Plant Biology, Uppsala BioCenter, Swedish University of 
Agricultural Sciences and Linnean Center for Plant Biology, SE-750 07 
Uppsala, Sweden. 2 Present address: Institute for Plant Protection, Italian 
National Research Council CNR, 50019 Sesto fiorentino, Italy, department of 
Crop Production Ecology, Swedish University of Agricultural Sciences and 
Linnean Center for Plant Biology, SE-750 07 Uppsala, Sweden, department 
of Plant Ecology and Evolution, Evolutionary Biology Centre, Uppsala 
University, SE-752 36 Uppsala, Sweden. 

Received: 27 September 2013 Accepted: 3 January 2014 
Published: 17 January 2014 



References 

1. Morgenstern EK: Geographic Variation in Forest Trees: Genetic Basis and 
Application of Knowledge in Silviculture. Vancouver, Canada: UBC Press; 1996. 

2. Levitt J: Growth and survival of plants at extreme temperature - A 
unified concept. In Dormancy and Survival, Symposia of the Society for 
Experimental Biology XXIII. Edited by Woolhouse HW. Cambridge: Cambridge 
University Press; 1969:395-448. 

3. Chuine I: Why does phenology drive species distribution? Phil Trans R Soc 
6 2010, 365:3149-3160. 

4. Chuine I, Beaubien EG: Phenology is a major determinant of tree species 
range. Ecol Lett 2001, 4:500-510. 

5. Nelson ND, Isebrands JG: Late-season photosynthesis and photosynthate 
distribution in an intensively-cultured Populus nigra x laurifolia clone. 
Photosynthetica 1983, 17:537-549. 

6. Weih M: Genetic and environmental variation in spring and autumn 
phenology of biomass willows {Solix spp.): effects on shoot growth and 
nitrogen economy. Tree Physiol 2009, 29:1479-1490. 

7. Pallardy SG: Physiology of Woody Plants. London, UK: Elseiver Inc.; 2008. 

8. Fuchigami LH, Weise CJ, Evert DR: Induction of cold acclimation in Cornus 
stolonifera Michx. Plant Physiol 1971, 47:98-103. 

9. Scarascia-Mugnozza GE, Hinckley TM, Stettler RF, Heilman PE, Isebrands JG: 
Production physiology and morphology of Populus species and their 
hybrids grown under short rotation. III. Seasonal carbon allocation 
patterns from branches. Can J Forest Res 1999, 29:1419-1432. 

10. Dodd RS, Huberli D, Mayer W, Harnik TY, Afzal-Rafli Z, Garbelotto M: 
Evidence for the role of synchronicity between host phenology and 
pathogen activity in the distribution of sudden oak death canker 
disease. New Phytol 2008, 179:505-514. 

11. Ghelardini L, Santini A: Avoidance by early flushing: a new perspective on 
Dutch elm disease. iForest 2009, 2:143-153. 

12. McKinney LV, Nielsen LR, Hansen JK, Kjaer ED: Presence of natural genetic 
resistance in Fraxinus excelsior (Oleraceae) to Chalara fraxinea 
(Ascomycota): an emerging infectious disease. Heredity 201 1, 106:788-797. 

13. Sinkkonen A, Somerkoski E, Paaso U, Holopainen JK, Rousi M, Mikola J: 
Genotypic variation in yellow autumn leaf colours explains aphid load in 
silver birch. New Phytol 2012, 195:461-469. 



14. Lagercrantz U: At the end of the day: a common molecular mechanism 
for photoperiod responses in plants? J Exp Bot 2009, 60:2501-2515. 

15. Howe GT, Gardner G, Hackett WP, Furnier GR: Phytochrome control of 
short-day-induced bud set in black cottonwood. Physiol Plant 1996, 
97:95-103. 

16. Bradshaw HD, Stettler RF: Molecular genetics of growth and development 
in Populus. IV. Mapping QTLs with large effects on growth, form, and 
phenology traits in a forest tree. Genetics 1995, 139:963-973. 

1 7. Keller SR, Soolanayakanahally RY, Guy RD, Silim SN, Olson MS, Tiffin P: Climate- 
driven local adaptation of ecophysiology and phenology in balsam poplar, 
Populus balsamifera L. (Salicaceae). Am J Bot 201 1, 98:99-108. 

18. Rohde A, Storme V, Jorge V, Gaudet M, Vitacolonna N, Fabbrini F, RuttinkT, 
Zaina G, Marron N, Dillen S, Steenackers M, Sabatti M, Morgante M, Boerjan 
W, Bastien C: Bud set in poplar - genetic dissection of a complex trait in 
natural and hybrid populations. New Phytol 201 1, 189:106-121. 

19. Weih M: Intensive short rotation forestry in boreal climates: present and 
future perspectives. Can J Forest Res 2004, 34:1369-1378. 

20. Tanino KK, Kalcsits L, Silim S, Kendall E, Gray GR: Temperature-driven 
plasticity in growth cessation and dormancy development in deciduous 
woody plants: a working hypothesis suggesting how molecular and 
cellular function is affected by temperature during dormancy induction. 
Plant Mol Biol 2010, 73:49-65. 

21. Rohde A, Bastien C, Boerjan W: Temperature signals contribute to the 
timing of photoperiodic growth cessation and bud set in poplar. 
Tree Physiol 2011,31:472-482. 

22. Heide OM, Prestud AK: Low temperature, but not photoperiod, controls 
growth cessation and dormancy induction and relaese in apple and 
pear. Tree Physiol 2005, 25:109-1 14. 

23. Heide OM: Interaction of photoperiod and temperature in the control of 
growth and dormancy of Prunus species. Sci Hortic-Amsterdam 2008, 
115:309-314. 

24. Heide OM: Temperature rather than photoperiod controls growth 
cessation and dormancy in Sorbus species. J Exp Bot 201 1, 62:5397-5404. 

25. Kozlowski TT, Pallardy SG: Growth Control in Woody Plants. San Diego, USA: 
Academic Press Inc; 1997. 

26. Fracheboud Y, Luquez V, Bjorken L, Sjodin A, Tuominen H, Jansson S: The 
control of autumn senescence in European aspen. Plant Physiol 2009, 
149:1982-1991. 

27. Keskitalo J, Bergquist G, Gardestrom P, Jansson S: A cellular timetable of 
autumn senescence. Plant Physiol 2005, 139:1635-1648. 

28. van der Schoot C, Rinne PLH: Dormancy cycling at the shoot apical 
meristem: transitioning between self-organization and self-arrest. 
Plant Sci 2011, 180:120-131. 

29. Perry TO: Dormancy of trees in winter. Science 1 971 , 171 :29-36. 

30. Polgar CA, Primack RB: Leaf-out phenology of temperate woody plants: 
from trees to ecosystems. New Phytol 201 1, 191:926-941. 

31. Heide OM: Dormancy release in beech buds (Fagus sylvatica) requires 
both chilling and long days. Physiol Plant 1993, 89:187-191. 

32. Sanz-Perez V, Castro-Diez P, Valladares F: Differential and interactive effects 
of temperature and photoperiod on budburst and carbon reserves in 
two co-occurring Mediterranean oaks. Plant Biol 2009, 1 1:142-151. 

33. Boucher LD, Manchester SR, Judd WS: An extinct genus of salicaceae based 
on twigs with attached flowers fruits, and foliage from the eocene green 
river formation of Utah and Colorado, USA. Am J Bot 2003, 90:1 389-1 399. 

34. Manchester SR, Judd WS, Handley B: Foliage and fruits of early poplars 
(Salicaceae : Populus) from the eocene of Utah, Colorado, and Wyoming. 
Int J Plant Sci 2006, 1 67:897-908. 

35. Berlin S, Lagercrantz U, von Arnold S, Ost T, Ronnberg-Wastljung AC: 
High-density linkage mapping and evolution of paralogs and orthologs 
in Salix and Populus. BMC Genomics 2010, 1 1:129. 

36. Kuzovkina Y, Weih M, Abalos Romero M, Belyaeva I, Charles J, Hurst S, 
Karp A, Labrecque M, Mclvor I, Singh NB, Smart LB, Volk TA: Salix: botany 
and global horticulture. Hortic Revs 2008, 34:447-489. 

37. Berlin S, Fogelqvist J, Lascoux M, Lagercrantz U, Ronnberg-Wastljung AC: 
Polymorphism and divergence in two willow species, Salix viminalis L. 
and Salix schwerinii E. Wolf. G3 201 1, 1:387-400. 

38. Junttila O: Apical growth cessation and shoot tip abscission in Salix. 
Physiol Plant 1976, 38:278-286. 

39. Junttila O: Effect of photoperiod and temperature on apical growth 
cessation in two ecotypes of Salix and Betula. Physiol Plant 1980, 
48:347-352. 



Ghelardini et al. BMC Plant Biology 2014, 14:31 
http://www.biomedcentral.com/1471-2229/14/31 



Page 1 8 of 1 8 



40. Ronnberg-Wastljung AC, Gullberg U: Genetics of breeding characters with 
possible effects on biomass production in Salix viminalis (L). Theor Appl 
Genet 1999, 98:531-540. 

41. Lennartsson M, Ogren E: Screening for efficient cold hardening in a 
breeding population of Salix using near infrared reflectance 
spectroscopy. Ann Forest Sci 2004, 61 :449-454. 

42. Tsarouhas V, Gullberg U, Lagercrantz U: Mapping of quantitative trait loci 
controlling timing of bud flush in Salix. Hereditos 2003, 138:172-178. 

43. Tsarouhas V, Gullberg U, Lagercrantz U: Mapping of quantitative trait loci 
(QTLs) affecting autumn freezing resistance and phenology in Salix. 
Theor Appl Genet 2004, 1 08:1 335-1 342. 

44. Pokhilko A, Fernandez AP, Edwards KD, Southern MM, Halliday KJ, Millar AJ: 
The clock gene circuit in Arabidopsis includes a repressilator with 
additional feedback loops. Mol Syst Biol 2012, 8:574. 

45. Hsu CY, Adams JP, Kim H, No K, Ma C, Strauss SH, Drnevich J, Vandervelde L, 
Ellis JD, Rice BM, Wickett N, Gunter LE, Tuskan GA, Brunner AM, Page GP, 
Barakat A, Carlson JE, DePamphilis CW, Luthe DS, Yuceer C: FLOWERING 
LOCUS T duplication coordinates reproductive and vegetative growth in 
perennial poplar. Proc Natl Acad Sci USA 201 1, 108:10756-10761. 

46. Ibanez C, Kozarewa I, Johansson M, Ogren E, Rohde A, Eriksson ME: 
Circadian clock components regulate entry and affect exit of seasonal 
dormancy as well as winter hardiness in Populus trees. Plant Physiol 2010, 
153:1823-1833. 

47. Dunlap JM, Stettler RF: Genetic variation and productivity of Populus 
trichocarpa and its hybrids.lX. Phenology and Melampsora rust 
incidence of native black cottonwood clones from four river valleys in 
Washington. Forest Ecol Manag 1996, 87:233-256. 

48. Frewen BE, Chen TH, Howe GT, Davis J, Rohde A, Boerjan W, Bradshaw HD 
Jr: Quantitative trait loci and candidate gene mapping of bud set and 
bud flush in Populus. Genetics 2000, 154:837-845. 

49. Howe GT, Saruul P, Davis J, Chen THH: Quantitative genetics of bud 
phenology, frost damage, and winter survival in an F-2 family of hybrid 
poplars. Theor Appl Genet 2000, 101:632-642. 

50. Fabbrini F, Gaudet M, Bastien C, Zaina G, Harfouche A, Beritognolo I, Marron 
N, Morgante M, Scarascia-Mugnozza G, Sabatti M: Phenotypic plasticity, 
QTL mapping and genomic characterization of bud set in black poplar. 
BMC Plant Biol 2012, 12:47. 

51. Marron N, Storme V, Dillen SY, Bastien C, Ricciotti L, Salani F, Sabatti M, Rae 
AM, Ceulemans R, Boerjan W: Genomic regions involved in productivity of 
two interspecific poplar families in Europe. 2. Biomass production and 
its relationships with tree architecture and phenology. Tree Genet 
Genomes 2010, 6:533-554. 

52. Barros RS, Neill SJ: Shoot growth in willow {Salix viminalis) in relation to 
abscisic acid, plant water status and photoperiod. Physiol Plant 1987, 
70:708-712. 

53. Molmann JA, Asante DKA, Jensen JB, Krane MN, Ernstsen A, Junttila O, Olsen JE: 
Low night temperature and inhibition of gibberellin biosynthesis override 
phytochrome action and induce bud set and cold acclimation, but not 
dormancy in PHYA overexpressors and wild-type of hybrid aspen. Plant Cell 
Environ 2005, 28:1579-1588. 

54. Kalcsits LA, Silim S, Tanino K: Warm temperature accelerates short 
photoperiod-induced growth cessation and dormancy induction in 
hybrid poplar (Populus x spp.). Trees 2009, 23:971-979. 

55. Chen THH, Howe GT, Bradshaw HD: Molecular genetic analysis of 
dormancy-related traits in poplars. Weed Sci 2002, 50:232-240. 

56. Bohlenius H, Huang T, Charbonnel-Campaa L, Brunner AM, Jansson S, 
Strauss SH, Nilsson O: CO/FT regulatory module controls timing of flowering 
and seasonal growth cessation in trees. Science 2006, 312:1040-1043. 

57. Ruttink T, Arend M, Morreel K, Storme V, Rombauts S, Fromm J, Bhalerao RP, 
Boerjan W, Rohde A: A molecular timetable for apical bud formation and 
dormancy induction in poplar. Plant Cell 2007, 19:2370-2390. 

58. Horvath DP, Chao WS, Suttle JC, Thimmapuram J, Anderson JV: 
Transcriptome analysis identifies novel responses and potential 
regulatory genes involved in seasonal dormancy transitions of leafy 
spurge (Euphorbia esula L). BMC Genomics 2008, 9:536. 

59. Yamane H, Ooka T, Jotatsu H, Hosaka Y, Sasaki R, Tao R: Expressional 
regulation of PpDAM5 and PpDAM6, peach (Prunus persica) 
dormancy-associated MADS-box genes, by low temperature and 
dormancy-breaking reagent treatment. J Exp Bot 201 1, 62:3481-3488. 



60. 



62. 



63. 



64. 



65. 



66. 



67. 



69. 



70. 



Diaz-Riquelme J, Lijavetzky D, Marti nez-Za pater JM, Carmona MJ: Genome-wide 
analysis of MIKCC-type MADS box genes in grapevine. Plant Physiol 2009, 
149:354-369. 

Li Z, Reighard GL, Abbott AG, Bielenberg DG: Dormancy-associated MADS 
genes from the EVG locus of peach [Prunus persica (L.) Batsch] have 
distinct seasonal and photoperiodic expression patterns. J Exp Bot 2009, 
60:3521-3530. 

Danfors B, Ledin S, Rosenqvist H: Short-rotation Willow Coppice. Growers' 
Manual. Uppsala, Sweden: Swedish Institute of Agricultural and 
Environmental Engineering; 1997. 

Van Ooijen JW: MapQTL ® 6, Software for the Mapping of Quantitative Trait Loci 
in Experimental Populations. Wageningen, Netherlands: Kyazama BY.; 2009. 
Lander ES, Botstein D: Mapping mendelian factors underlying 
quantitative traits using RFLP linkage maps. Genetics 1989, 121:185-199. 
Churchill GA, Doerge RW: Empirical threshold values for quantitative trait 
mapping. Genetics 1994, 138:963-971. 

Lander E, Kruglyak L: Genetic dissection of complex traits: guidelines for 
interpreting and reporting linkage results. Nat Genet 1995, 11:241-247. 
Van Ooijen JW: LOD significance thresholds for QTL analysis in 
experimental populations of diploid species. Heredity 1999, 83:613-624. 
Jansen RC, Stam P: High resolution of quantitative traits into multiple loci 
via interval mapping. Genetics 1994, 136:1447-1455. 
Snijders TAB, Bosker RJ: Multilevel Analysis: An Introduction to Basic and 
Advanced Multilevel Modeling. London, UK: Sage Publications Ltd; 1999. 
R Development Core Team: R: A Language and Environment for Statistical 
Computing. Vienna, Austria: R Foundation for Statistical Computing; 2010. 



doi:1 0.1 1 86/1 471 -2229-1 4-3 1 

Cite this article as: Ghelardini et al.: Genetic architecture of spring and 
autumn phenology in Salix. BMC Plant Biology 2014 14:31. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



